Metamaterials with conformational nonlinearity 
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Within a decade of fruitful developments, metamaterials became a prominent area of research, 
bridging theoretical and applied electrodynamics, electrical engineering and material science. Being 
man-made structures, metamaterials offer a particularly useful playground to develop novel interdis- 
ciplinary concepts. Here we demonstrate a novel principle in metamaterial assembly which integrates 
electromagnetic, mechanical, and thermal responses within their elements. Through these mecha- 
nisms, the conformation of the meta-molecules changes, providing a dual mechanism for nonlinearity 
and offering nonlinear chirality. Our proposal opens a wide road towards further developments of 
nonlinear metamaterials and photonic structures, adding extra flexibility to their design and control. 



Metamaterials — artificial materials designed to de- 
liver an unusual electromagnetic response — have already 
established a prominent area of theoretical and experi- 
mental physics 1 2 with applications ranging from super- 
imaging^ and transformation optics^ to tunable and ac- 
tive materials 5 , photonics^ and plasmonics^. At the same 
time, the corresponding ideas were readily adopted in 
acoustics, promising unusual mechanical behavioural 
Clearly, simultaneous access to electromagnetic, me- 
chanical and thermal properties offers a great applica- 
tion capabilities, as demonstrated, e.g. by mechanically 
tunable^ and thermally reconfigurable 11 metamaterials. 
Such a connection has been so far achieved by engineer- 
ing the metamaterial structure so as to provide an ar- 
tificial means to control electromagnetic properties; this 
approach is useful for tunability, however, it does not 
provide a dynamic interaction. 

Here we design a simple metamaterial element where 
responses of a different nature are intrinsically coupled 
through its very structure. Such an element is readily 
found as a thin wire wound into a spiral (Fig. [T]). From 
an electromagnetic perspective, it is a resonator with chi- 
ral properties, and from a mechanical perspective it is a 
spring. Finally, its specific shape makes all the character- 
istics sensitively depending on temperature. This natural 
"coincidence" is of a great use as soon as the parameters 
of the spiral allow thermo-magneto-mechanical coupling, 
which results in efficient changes of its conformation. In- 
deed, an external electromagnetic wave induces a current 
along the spiral conductor, but this current causes an 
attractive force between the windings. The spiral will 
therefore contract until that force is balanced by the 
spring force, however the corresponding change in spi- 
ral geometry will shift the resonance and thus alter the 
current amplitude in a self-consistent manner. In addi- 
tion, the thermal expansion of the spiral, resulting from 
its heating with the increasing incident power, provides 
a relevant contribution, further shifting the resonance — 
importantly, in the same direction. This dual conforma- 
tional feedback leads to exotic nonlinear behaviour. 

Spirals are chiral resonant particles known for about a 
centur^l, being suitable for producing bi-isotropid^ or 
bi-anisotropicP^ media, depending on how they are ar- 
ranged in space. Approximate models were developed 
to derive polarizability of simplified canonical spirals^, 
which can be used to derive effective parameters of di- 
lute media, as well as for bi-anisotropic lattices of infi- 
nite spirals 16 . However, in the general case, the analy- 
sis of a periodic or random lattice of volumetric spiral 




FIG. 1. Schematic of a metamaterial composed of spiral 
meta-molecules, which are electromagnetic resonators which 
can change their geometrical conformation: spiral pitch £ can 
vary, following a balance between the attractive force induced 
by magnetic field and the spring rigidity; and spiral radius r 
can change upon thermal expansion. 



scatterers eludes analytical treatment. For this reason, 
in the qualitative analysis in this paper we assume that 
the individual elements are arranged in such a way as to 
minimize mutual interaction so that the mutual coupling 
would only manifest itself in a small shift of the resonance 
frequency^. We can therefore retrieve all the qualitative 
features by analysing the response of an individual spiral. 

Nonlinear behaviour of the metamaterials made of spi- 
ral meta-molecules can be theoretically understood as 
its resonance frequency essentially depends on the spiral 
pitch £ as well as on the radius of the spiral r. Accord- 
ingly, the resonance shifts when the spiral compresses in 
response to attraction between the windings as well as 
when a thermal expansion occurs due to heating, both 
the effects triggered by the induced currents. This en- 
sures the intensity-dependent phenomena, providing for 
essential nonlinearity, for example, in dependence of the 
induced magnetization on the incident magnetic field. 
We note that the change in spiral pitch £ manifests itself 
with the corresponding change of its chiral properties. 
We can qualitatively describe the degree of chirality with 



the ratio between the electric p z and magnetic m z mo- 
ments, induced in each spiral along its axis z. This ratio 
is directly determined by the geometry of the spiraP^, 
\l\ \p z /m z \ f/Wr. 

Precise electromagnetic modelling of a multi-turn finite 
spiral is known to be analytically challenging, however for 
a good qualitative description of the nonlinear phenom- 
ena, we can rely on a simple circuit modet^ which is 
applicable to two-windings spiral resonators if their res- 
onance frequency is well within the quasi-static regime 
(see Methods). 

For small £, the Ampere force acting between the wind- 
ings of the spiral, can be calculated as that between the 
two parallel wires of the corresponding length. Given the 
actual current distribution along the entire spiral, we 
have found the net total force acting between the wind- 
ings to be F u = This attractive force is be 
countered by the spring force F s (see Methods). The me- 
chanical equilibrium is thus governed by 



3GV 2 «, 4 (£ - + 2/zo/ 2 = 0, 



(1) 



which may seem to be a quadratic equation for £, however 
it is in fact far more complicated as the current I also 
depends on £ and r through the impedance equation: 



R(r) + iuL(r) 



• / = —\ujfiQTir 2 HQ. (2) 



Here, Ho is the amplitude of the magnetic field of the 
incident wave (we imply the incident polarization with 
H being parallel to the spiral axis). 

The thermal behaviour of the spiral is provided by its 
thermal expansion r(T) as well as by the temperature de- 
pendence of its conductivity (see Methods). The amount 
of heat Q + = RI 2 /2 which is realized per second in the 
spiral, is balanced, in a thermal equilibrium, to the heat 
loss through its surface Q~ = j3 • 2w(2i\r) 2 • AT, where /3 
is a an empirical coefficient which depends on the cooling 
conditions. Thus, equations ([I]), Q and 



P ■ 2w(2irr) 2 ■ AT = RI 2 /2 



(3) 



form a system of coupled equations, which can be numer- 
ically solved for £, AT and I at a given frequency uo and 
amplitude Ho of the incident field. 

For a physical representation of the arising nonlinear- 
ity, we can evaluate the total magnetic moment, induced 
in each spiral, and obtain the effective magnetization of 
the metamaterial, M = Ivuro (where v is the volumetric 
concentration of spirals). 

In Fig. [2]a-c, we illustrate the amplitude-dependent 
phenomena observed when the signal frequency is slightly 
lower than uoq. With the increase of power, the spiral 
compression enhances the resonance, and therefore a pos- 
itive feedback occurs upon a certain threshold, when the 
spring force becomes insufficient to hold the resonantly 
growing attraction [follow the green right-pointing trian- 
gles in Fig. [2|i]. The spiral then jumps to another con- 
formation, where a stronger compression shifts the reso- 
nance sufficiently below the signal frequency, so that the 
induced force becomes weak again, and further increase 
in power causes slow growth in magnetization. However, 
if the field intensity is decreased from the high values, the 
spiral remains stable in a compressed conformation much 
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FIG. 2. Dependence of the metamaterial properties on either 
(a-c) the incident amplitude Ho (at the relative frequency of 
0.99cjo), or (d-f) on the signal frequency (with the amplitude 
2 A/m), for relaxed (right-pointing green triangles) and com- 
pressed (left-pointing red triangles) conformations. Panel (a) 
shows the change in spiral pitch A£; panel (b) the increase 
in spiral temperature AT, panel (c) the absolute value of the 
effective magnetization, panel (d) the relative shift in the res- 
onance frequency Acj, panel (e) the magnetic moment (real 
and imaginary parts) of a single spiral, and panel (f) the chi- 
rality parameter 7 • cq . 



longer, as it is now trapped near the resonance. It is only 
released when the peak amplitude of force is insufficient 
to hold it [follow the red left-pointing triangles in Fig. [2^1]. 
Further contribution is provided by the thermal expan- 
sion of the spiral heated by the currents (Fig. [2]o) . The 
two effects act simultaneously, providing a remarkable 
hysteresis-like pattern for all the macroscopic properties 
of the metamaterial, such as magnetization (Fig. [2^). 

Two conformations can be also observed when the fre- 
quency of the incident radiation is varied (Fig. |2]d-f). 
When the signal frequency increases at a given ampli- 
tude, the spiral remains just slightly compressed until, 
close to the initial resonance c^o, the attractive force be- 
comes strong enough to trigger a jump to a compressed 
conformation (so that the spiral resonance goes below the 
signal frequency), and further frequency increase slowly 
relaxes this compression. With decreasing frequency, 
the resonance of the spiral is already below the sig- 
nal, so approaching it steadily increases attractive forces 
and the spiral undergoes continuous compression, led by 
the slope of the resonance. To illustrate the frequency- 
dependent phenomena, we show how the resonance fre- 
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FIG. 3. Experimental results. Panel (a) shows transmis- 
sion through the spiral as a function of frequency and input 
power, panel (b) the curves for the highest and lowest mea- 
sured power, panel (c) the frequency shift as a function of 
input power. 
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quency (Fig. the magnetic moment (Fig. [2p 
the chirality (Fig. [2]p) vary in hysteresis-like loops. 

For the proof-of-principle experiments in the mi- 
crowave range, we take a thicker wire, such that the ther- 
mal effect dominates. Fig. [3](a) shows the transmission 
amplitude as a function of input frequency and power. 
It can be seen that there is a substantial shift of the 
resonance to lower frequencies as the input power is in- 
creased. As can be seen in Fig. [3jb) , the maximum fre- 
quency shift is less than the width of the resonances, 
which is insufficiently strong to observe bistability. Fig- 
ure [3|c) shows a shift of the resonant frequency as a func- 
tion of the input power. Analysis of this shift allows to 
estimate the coefficient in Eq. |3| which determines the 
effective heat loss from the spiral, which turns out to be 
about f3 « 30. We have used this coefficient for the above 
theoretical calculations. 

The above results demonstrate that the flexible spiral 
resonator exhibits nonlinear behaviour with two stable 
conformations of its meta-molecules. Clearly, to ensure 
sufficiently strong force between the current loops, the 
windings of the spiral have to be close to each other, 
£ <C 1, while the wire should be thin enough to make 
the spring soft to respond to the relatively weak Ampere 
forces. 

Among the parameters of the spiral, the relative wire 
thickness w has the most drastic effect on the interplay 
between the mechanical response (which increases as the 
fourth power of 1/w) and the thermal effect (largely inde- 
pendent of w). We conclude that for the spiral's compres- 
sion to dominate over its thermal expansion the relative 
wire width should be smaller than 0.01. The conforma- 
tional dynamics of the spiral is most interesting when the 
thermal and mechanical reactions are comparable to each 
other, as shown with our theoretical examples. 

Our proposal provides a starting platform for further 
analysis of the nonlinear phenomena in various metama- 



terials based on the flexible chiral particle, which can be 
made bi-isotropic, bi-anisotropic and even non-chiral, by 
choosing an appropriate lattice and alignment. While a 
lot of peculiar features will be provided with the arising 
mutual interactions, the essence of nonlinear response is 
determined by the properties of an individual spiral par- 
ticle, and therefore the above analysis forms a reliable 
basis for further research as well as the development of 
useful applications. 



METHODS 

In our theoretical analysis, r is the radius of the spiral loops, 
r w is the radius of the wire and d is the vertical distance between 
windings; the dimensionless normalized quantities w = r w /r and 
£ = d/r are used for convenience (Fig. 

With the quasi-static models the spiral is described as an LC- 
circuit, with the inductance L being the same as the inductance of a 
single circular loop of the same wire, L = fior (ln(8 / w) — 2) (because 
the net current, flowing within the entire spiral — the sum over 
the current distribution in all turns — is uniform along the entire 
circumference). The latter equally applies to the resistance R = 
a/ noou/2a/w] both the L and R are taken with skin-effect in mind. 
The capacitance C is defined by that between the two windings 
(with a minor correction for their curvature being negligible for 
small £): C(£) = 2tt£q • 7rr/acosh(£/2K;). We have checked, with 
the aid of numerical simulations, that this simple circuit model 
correctly predicts the resonant frequency uv(£) = 1/ \/L(r)C(r, £) 
of two-turn spirals (Fig.^J. 

The mechanical properties of a two- windings spiral are described 
by the stiffness coefficient k = Griu 4 /8, where G is the shear modu- 
lus of the spiral material. Note that the characteristic frequency of 
mechanical oscillations, uom — w^3G/2p/ (47rr) is many orders of 
magnitude smaller than the electromagnetic frequencies involved, 
so the analysis of the spiral geometrical reconfiguration is essen- 
tially static and is determined by time-averaged amplitudes of the 
current. The spring response is then described with the Hooke's 
law, so it linearly increases with the deviation from the initial pitch 
value £ : = kr(^ — fo)- 

Temperature dependence of the spiral parameters is determined 
by a thermal expansion coefficient a as r = ro(l-\-a-AT). Note that 
it leaves £ and w unchanged provided that the spiral is perfectly 
annealed and is free from internal stress. The electric conductivity 
is affected in a similar way, a = (Jq(1 + 77 • AT) with a temperature 
coefficient of conductivity 77. 

For theoretical calculations, we assumed rrj = 5 mm, £0 = 0.04, 
w = 0.01, and the material parameters of copper for a = 6.6- 10 7 S, 
G = 40GPa, a = 1.7- lO" 5 K" 1 and 77 = 4.3 • lO" 3 K _1 . 

Numerical simulations were performed using the finite-element 
based frequency-domain solver of CST Microwave Studio. The spi- 
ral was excited by a plane wave, with its magnetic field parallel to 
the axis of the the spiral. The material was modelled as a perfect 
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FIG. 4. Resonance frequency of spiral resonators with differ- 
ent pitch £. Comparison between the effective circuit theory 
(solid line) and numerical simulations (circles). Spiral param- 
eters here are ro — 2 mm and w — 0.01. 
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electrical conductor, and a curved mesh was used to ensure good 
resolution of the geometry without excessively fine meshing. The 
spiral was excited over a broad frequency range, and the resonant 
frequency was determined as that with the maximum magnitude 
of the normal magnetic field at the centre of the spiral. 

For our experiments, the spirals were made of 100 /im-thick cop- 
per wire with 2.4 mm radius and 0.5 mm pitch. We suspend a 
spiral parallel to the axis of a WR 229 rectangular waveguide on a 
dielectric rod so that the incident magnetic field is along the spi- 
ral axis. The waveguide is excited by a vector network analyzer 
(Rhode and Schwartz ZVB-20) after amplification by a 1W am- 
plifier (HP 83020A). The output of the waveguide is measured by 
the network analyzer, through a 20 dB attenuator to eliminate gain 
compression in the detectors. To ensure that the nonlinear response 
of the amplifier does not contribute to the measurements, a power 
calibration is performed for each individual power level, thus en- 



suring that the power level at the waveguide input is compensated 
for the frequency and power dependent gain of the amplifier. A 
directional coupler is used to sample the output from the amplifier, 
so that the power calibration can be performed in situ, and to pro- 
vide a reference level against which to compare the output signal 
from the waveguide. The input power to the waveguide was varied 
between 5dBm and 28dBm in steps of ldB. 
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